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Abstract 



^SJ ■ We show that an ra-gap periodic quantum system with parity-even smooth potential 

admits 2 n — 1 isospectral super-extensions. Each is described by a tri-supersymmetry 
that originates from a higher-order differential operator of the Lax pair and two-term 
nonsingular decompositions of it; its local part corresponds to a spontaneously par- 
tially broken centrally extended nonlinear N = 4 supersymmetry. We conjecture that 
any finite-gap system having antiperiodic singlet states admits a self-isospectral tri- 
supersymmetric extension with the partner potential to be the original one translated 
CT} ■ for a half-period. Applying the theory to a broad class of finite-gap elliptic systems de- 

scribed by a two-parametric associated Lame equation, our conjecture is supported by 
the explicit construction of the self-isospectral tri-supersymmetric pairs. We find that 
the spontaneously broken tri-supersymmetry of the self-isospectral periodic system is 
recovered in the infinite period limit. 

O 

1 Introduction 



Finite-gap periodic quantum systems play an important role in physics. They underly the 
theory of periodic solutions in nonlinear integrable systems, including the Korteweg-de Vries, 
the nonlinear Schrodinger, the Kadomtsev-Petviashvili, and the sine-Gordon equations [1, 
2, 4, 3, 5]. Being analytically solvable systems, they find various applications in diverse 
areas. The list of their applications is extensive, and among others includes the modelling 
of crystals [6, 7, 8], the theory of monopoles [9], instantons and sphalerons [10, 11], classical 
Ginzburg-Landau field theory [12], Josephson junctions theory [13], magnetostatic problems 
[14], inhomogeneous cosmologies [15], Kaluza- Klein theories [16], chaos [17], preheating 
after inflation modern theories [18], string theory [19], matrix models [20], supersymmetric 
Yang-Mills theory [21, 22] and AdS/CFT duality [23]. 

Some time ago it was showed by Braden and Macfarlane [24] , and in a more broad context 
by Dunne and Feinberg [25], that a usual N = 2 supersymmetric extension of a periodic 
quantum system may produce a completely isospectral pair with a zero energy doublet of 
the ground states. Such a picture is completely different from that taking place in non- 
periodic systems described by the same linear N = 2 super algebraic structure {Q a , Qb} = 
25 a bH, [Q a , H] = 0. There, the complete isospectrality of the super-partners happens only 
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in the case of a spontaneously broken supersymmetry, characterized by a positive energy 
of the lowest supersymmetric doublet [26]. Furthermore, it was showed that there exist 
peculiar isospectral supersymmetric periodic systems, in which the partner potentials are 
identical in shape but mutually translated for a half-period, or reflected, or translated and 
reflected. A pair of super-partner potentials with such a property was named by Dunne 
and Feinberg as self-iso •spectral. The phenomenon of self-isospectrality with a half-period 
shift was illustrated by some examples of exactly soluble models belonging to a class of 
finite-gap periodic systems. Later on isospectral and self-isospectral supersymmetric finite- 
gap periodic systems were studied in various aspects [27, 28, 29], and it was found in [30] 
that a property of the self-isospectrality may also appear in some periodic finite-gap systems 
based on a nonlinear supersymmetry of the second order {Q a , Q b } = 2<5 afe P 2 (-£0, [Q a , H] = 
[31, 32, 33], where P2(H) is a quadratic polynomial. The nature and origin of isospectrality 
and self-isospectrality in finite-gap systems have remained, however, to be obscure. 

Recently we showed [34] that self-isospectrality may be realized by a non-relativistic 
electron in the periodic magnetic and electric fields of a special form, and indicated on a 
peculiar non-linear supersymmetric structure associated with it. 

In the present paper, we study superextension of quantum periodic systems with a parity- 
even smooth finite-gap potential of general form, and show that it is characterized by an 
unusual tri-supersymmetric structure. This peculiar supersymmetric structure originates 
from the higher order differential operator of the Lax pair, and its decomposability in pairs 
of nonsingular operators. The superalgebra, generated by three indicated integrals of motion 
together with trivial integrals associated with parity symmetry and matrix extension, has a 
nonlinear nature, that reflects a nonlinearity of a spectral polynomial of the original finite-gap 
system. The higher order operator of the Lax pair of a nontrivial n-gap (n > 0) system admits 
2 n — 1 nonsingular two-term decompositions. By means of the Crum-Darboux transformation, 
with each nonsingular decomposition we associate a particular tri-supersymmetric extension, 
and as a result get a family of 2 n completely isospectral systems. We show that a local part 
of the tri-supersymmetry is a spontaneously partially broken centrally extended nonlinear 
N = 4 supersymmetry, that explains the nature and origin of the complete isospectrality. 
When the original finite-gap system has in its spectrum a nonzero number of anti-periodic 
singlet states corresponding to the edges of permitted bands, among all the non-singular 
decompositions of the higher order operator of the Lax pair there is a special one which 
corresponds to a separation of all the singlets into orthogonal subspaces of periodic and anti- 
periodic states. We conjecture that it is this separation that produces a self-isospectral tri- 
supersymmetric system. This means particularly that all the set of 2 n completely isospectral 
systems we get, including the original n-gap system with the specified special property, is 
divided into 2 n_1 self-isospectral tri-supersymmetric pairs. Then we apply a general theory 
to a broad class of finite-gap elliptic (double periodic) quantum systems described by a two- 
parametric family of associated Lame equations. Any such a system has in its spectrum a 
non-empty subspace of anti-periodic singlet states, and we support our conjecture by the 
explicit construction of the self-isospectral tri-supersymmetric pairs. We also investigate a 
rather intricate picture of the infinite-period limit of the tri-supersymmetry. 

The paper is organized as follows. In the next section we first discuss general properties 
of the finite-gap periodic systems with a smooth potential, and show that any parity-even n- 
gap system is characterized by a hidden bosonized N = 2 nonlinear supersymmetry of order 
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2n + 1. This supersymmetry reflects the peculiarities of the band structure. In Section 3 we 
show how the tri-supersymmetric extensions of the system are constructed by means of the 
Crum-Darboux transformation. There we also investigate a general structure and properties 
of the tri-supersymmetry. In Section 4 we apply a general theory to the case of finite-gap 
elliptic systems described by the associated Lame equation. In section 5, the infinite-period 
limit of the tri-supersymmetry is studied. Section 6 is devoted to concluding remarks and 
outlook. 



2 Hidden supersymmetry in finite-gap systems 

To have a self-contained presentation, in this section we first summarize briefly the properties 
of the quantum periodic systems of a general form. Then we restrict the consideration to 
the case of the smooth parity-even finite-gap systems to reveal in them a hidden bosonized 
nonlinear supersymmetry whose order is defined by the number of energy gaps. 

2.1 General properties of quantum periodic systems 

Consider a quantum system given by a Hamiltonian operator H = —D 2 + u(x), D = with 
a real smooth periodic potential u(x), u(x) = u(x + 2L). For the corresponding stationary 
Schrodinger equation, 

HV(x) = EV(x), (2.1) 

known in the literature as Hill's equation, we choose some real basis of solutions, ipi(x; E), 
ip2{x] E). The operator of translation for the period 2L, or the monodromy operator, 

TV(x) = ^(x + 2L), (2.2) 

commutes with the Hamiltonian H, [T, H] = 0. It preserves a two-dimensional linear vector 
space of solutions of (2.1), and can be represented there by the second order monodromy 
matrix M(E), 

T^ a (x; E) = ^ a (x + 2L; E) = M ab (E)ip h (x; E). (2.3) 

The change of the basis, t/j a (x; E) — > ip a (x; E) — A ab ^ b (x; E), det A ^ 0, generates a con- 
jugation of the monodromy matrix, M(E) — > M(E) = AM(E)A~ 1 , but do not change its 
determinant, det M(E) = det M(E), trace, TrM(E) = TrM = V(E), and eigenvalues, given 
by solutions of the characteristic equation 

det(M(£) - fil) = 0. (2.4) 

Let us choose a particular basis of solutions fixed by conditions 

M0;E) = 1, ,f>[(p;E)=0, M0;E)=0, ^' 2 (0;E) = 1, (2.5) 

where prime denotes the x-derivative. Differentiating relation (2.3) in x and putting then 
x = in (2.3) and in the derived relation, we find the form of the monodromy matrix in 
basis (2.5), 

( ^(2L;E) ^(2L;E) \ 
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Wronskian W(ip\, ^2) = ^1^2 — ^1^2 of any two linearly independent solutions of equation 
(2.1) takes a nonzero ^-independent value, which for basis (2.5) is equal to 1. Then the 
explicit form of the real monodromy matrix (2.6) shows that a basis-independent value of 
its determinant does not depend on energy either, det M(E) = 1, and so, M(E) G sl(2,M). 
Note that the change x = — > rro G K in relations (2.5) gives a one-parametric family of 
the bases, i[> a (x; xo, E) = A(xo) a bipb(x; E), A{xq) G s/(2,M), playing an important role in 
the theory of periodic quantum systems [l]-[4]. In such a basis, the monodromy matrix will 
include an additional dependence on the marked point xq, M(E,xq) G s/(2,M). 

With taking into account that detM = 1, the characteristic equation (2.4) is reduced to 
1 — V{E)ji + /i 2 = 0, and the basis-independent eigenvalues of the monodromy matrix are 
given in terms of its trace 1 , 



In correspondence with det M(E) = 1, /ii/z 2 = 1- Common eigenstates of H and T are 
described by the Bloch-Floquet functions ip±(x; E), which satisfy a relation 



The values of the discriminant T>(E) define the spectral properties of the periodic 
Schrodinger equation. For some energies E G (-E^j-i, E 2 i), % — 0, . . ., Ei < E i+1 , E_i = — 00, 
the quasi-momentum k(E) takes complex values, and \V(E)\ > 2. Solutions corresponding 
to such £"s are not physically acceptable as they diverge in x = —00 or +00. For these 
values of E we have a forbidden band, or an energy gap, see Fig. 1. In a generic case, a 
periodic quantum system has an infinite number of gaps. The width of the gaps decreases 
rapidly when energy increases, while the rate of decrease depends on the smoothness of the 
potential. In the case of analytic potentials, the gaps decrease exponentially. Energies E for 
which |X>(i?)| < 2, define permitted bands, or permitted zones. Here, the quasi-momentum 
k(E) takes real values, and complex numbers exp(±iK,(E)) have modulus equals to 1. All 
the energy levels with \V(E)\ < 2 are doubly degenerate, but for \V(E)\ = 2 we have two 
essentially different cases. For those E, which separate permitted and prohibited bands, 
corresponding eigenvalue of the monodromy matrix is non-degenerate, the matrix M has a 
form of Jordan matrix, and a physical singlet band-edge state is periodic, exp(iK,(E)) = +1, 
if T>(E) = 2, while for T>(E) = —2 a singlet state is antiperiodic, exp(iK,(E)) = —1. When 
|r>(-E)| = 2 but the corresponding eigenvalue of the monodromy matrix is doubly degenerate, 
M is diagonalizable on the two linearly independent Bloch-Floquet states, which both are 
periodic if T>(E) = 2, or are antiperiodic when T>(E) = —2. This second situation, that 
corresponds to points E 3 = E4 and Eg = E w on Fig. 1, takes place when a prohibited band 
disappears. 

Summarizing, the interval (—00, E ) constitutes the lowest forbidden band. The permit- 
ted bands with \ V(E)\ < 2 are separated by prohibited bands, or energy gaps. All the energy 

lr The trace of the monodromy matrix is called in the literature the Lyapunov function, Hill determinant, 
or discriminant of the Schrodinger equation. 




(2.7) 



T^ ± (x; E) = exp(±iK(E))tp±(x] E), 
where /j,i t2 (E) = exp(±in(E)), and the quasi-momentum k(E) is given by 

2cos/<£) = V(E). 



(2.9) 



(2.8) 
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Figure 1: The discriminant V{E) in a generic situation of a periodic potential. 



levels in the interior of permitted bands have a double degeneration, while the states at their 
edges are singlets. 

According to the oscillation theorem [35], the common eingenstates of H and mon- 
odromy operator T with energies E < E 1 < E 2 < E 3 < E 4 < E 5 < E 6 < . . . such 
that \V(Ek)\ = 2, are described by the wave functions which are characterized by the pe- 
riods 2L, 4L, 4L, 2L, 2L, 4L, 4L..., and by the node numbers in the period 2L equal to 
0, 1, 1, 2, 2, 3, 3, . . ., see Fig. 1. The odd number of nodes corresponds to antiperiodic states, 
whereas the periodic states have an even number of nodes in the period 2L. The singlet 
states at the edges of the same prohibited band have the same number of nodes and the 
same periodicity, and their nodes are alternating. 

2.2 Finite-gap systems and hidden bosonized supersymmetry 

In some periodic potentials infinite number of bands merge together so that only finite 
number of gaps remains in the spectrum. Such potentials are called finite-gap. The simplest 
case of a zero-gap system corresponds here to a free particle with u(x) = const 2 . For the 
Schrodinger equation with a finite-gap potential the spectrum and eigenf unctions can be 
presented in an analytical form 3 . Having also in mind that for analytical potentials the 
size of the gaps decreases exponentially when energy increases, any periodic potential can 
be approximated by a finite-gap potential if narrow gaps are disregarded. 

From now on we suppose that a periodic potential u(x) is finite-gap. Additionally, we 
assume that it is an even function, u(x) = u(—x). Then a reflection (parity) operator 
R, Rip{x) = ip(—x), is a nonlocal integral of motion, [R, H] = 0. Periodicity and parity 
symmetry together imply that the potential possesses also a middle-point reflection symmetry 
u(L + x) — u(L — x). 

2 Here and in what follows we do not count the prohibited band (— oo, E n ) that always presents in any 
periodic system with a smooth potential. 

3 In this sense, and in a contrast with, for example, the Kronig-Penney model, finite-gap potentials play 
the same role in solid-state physics as the Kepler problem in atomic theory. 
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The spectrum a(H) of a nontrivial n-gap (n > 0) system is characterized by the band 
structure, a(H) = [E , E x ] U . . . U [E 2n _ 2 , E 2n -i\ U [E 2n ,oo), where E < E 1 < . . . < E 2n 
are the non-degenerate energies corresponding to the 2n + 1 singlet band-edge states ^i(x), 
H^i = E^i, i = 0,1,..., 2n. Since parity operator R is an integral, each singlet state 
^i(x) has a definite parity, +1 or —1. The energy levels in the interior of permitted bands, 
E G (E 2i , E 2i+ i), % = 0, ...,n, are doubly degenerate, and certain linear combinations of 
corresponding Bloch-Floquet doublet states are the eigenstates of R with eigenvalues +1 and 
— 1. These properties indicate on the presence of a hidden bosonized N = 2 supersymmetry in 
any finite-gap system, for which operator R has to play the role of the grading operator. The 
presence of 2n + 1 > 3 singlet states indicates, however, on its nonlinear nature [31, 32, 33]. 
The supercharges and the form of the corresponding nonlinear superalgebra can easily be 
identified. 

Any finite-gap system is characterized by the presence of a nontrivial integral of motion 
in the form of an anti-Hermitian differential operator of order 2n + 1, 

A 2n+1 = D 2 " +1 + 4(x)D^ + 4(x)D 2n - 2 + ... + ct(x) : (2.10) 

where the coefficient functions cf(x) are real. The absence of the term proportional to D 2n 
in its structure, i.e. an equality c±(x) = 0, is dictated by the condition [A 2n+ i,H] = 0. 
Other coefficients cf(x) are fixed in the form of polynomials in the potential u(x) and its 
derivatives [5]. Thus, for periodic potential, A 2n+ i is a periodic operator, i.e. [A 2n+ i,T] = 0. 
(A 2n+ i, H) is known as the Lax pair of the n-th order Korteweg-de Vries (KdV) equation. 
A possible form of the n-gap potential is fixed by a nonlinear equation, which has a sense of 
the n-th equation of the stationary KdV hierarchy [2, 5]. This equation can be represented 
alternatively as 

L(JL)"1 = 0, L = D 3 + 2 (uD + Du), (2.11) 

where J is the operator of indefinite integration J = D~ l [4]. The form of a one-gap 
potential is fixed by this equation in a unique manner, u(x) = TP[x + oj 2 + c), where V(x) is 
the doubly periodic (elliptic) Weierstrass function [36], and c is a constant. To have a real- 
valued potential, one of the periods of V(x) is chosen to be real, 2ui = 2L, while another 
period 2u 2 is assumed to be pure imaginary, and cGl In the case n > 1 the form of the 
potential u(x) is not fixed uniquely even if it is restricted to a class of elliptic functions. 
The mutually commuting operators A 2n+ i and H satisfy the relation 

2n 

-A 2 2n+1 = P 2n+1 (H), P 2n+1 (H) = ]J(H - Ej), (2.12) 

3=0 

where P 2n+ i(H) is a spectral polynomial given in terms of singlet energies. It is in accordance 
with Burchnall-Chaundy theorem [37, 38], which says that if two differential in x operators 
A and B of mutually prime orders / and k, respectively, commute, [A, B] = 0, they satisfy a 
relation P(A, B) = 0, where P is a polynomial of order k in A, and of order / in B. Equation 
(2.12) corresponds to a non-degenerate (Ei ^ Ej for i ^ j) spectral elliptic curve of genus n 
associated with an n-gap periodic system [l]-[4] 4 . 

4 Because of the described properties, u(x) is called algebro-geometric finite-gap potential. 
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As a consequence of (2.12), the operator A 2n +i annihilates all the 2n + l singlet band-edge 
states. Indeed, from [A 2n +i, H] = we have A 2n+ i^j = a^j + P&j, where is a physical 
{T^/j = 'y'&j, 7 G { — 1, 1}) and $j is a non-physical solution corresponding to a band-edge 
energy Ej. Acting from the left by T, we get r yA 2n +i^j = 'ja^/j + (3TQj, and, therefore, 
(3(jT — = 0. As $j is neither periodic nor antiperiodic, the last equation can be satisfied 
if and only if f3 — 0. Then, equation (2.12) dictates that a — 0. 

Consider the Wronskian of the singlet states, W A = W^(^ , ^2n)- In a generic case 
the Wronskian of s linearly independent functions that form a kernel of an arbitrary linear 
differential operator of order s, C = D s + ci(x)D s ^ 1 + . . ., satisfies the Abel identity W'(x) = 
— ci(x)W [38]. For operator (2.10) a corresponding coefficient function is c A {x) = 0, and 
because of the linear independence of the singlet band-edge states we find that 

W A (x) = C + 0, (2.13) 

where C is a constant. When s linearly independent zero modes ipj, j — 1, . . . , s, of operator 
C are known, the form of this operator can be reconstructed in their terms. The coefficients 
Ck(x) are defined by relations Ck(x) — — ^r, k — 1, . . . , s, where the functions Wk(x) are 

obtained from Wronskian W = W((pi, . . . ,tp s ) by replacing in it y^ 5 ^ = D s ~ h (fj by (p^ 
[40] , see Appendix A. In our case, each singlet band-edge state ^i(x), being a zero mode of 
A 2n +i, possesses a definite parity. As a result, with taking into account (2.13), we find that 
the coefficients c^ r (x) are odd, while the coefficients c^ r+1 (x) are even non-singular functions. 
Hence the integral A 2n+ i is parity odd, 

{R,A 2n+1 } = 0. (2.14) 

Introducing two Hermitian operators 

Z = Z x = iA 2n+1 , Z 2 = iRZ, (2.15) 

and identifying them as odd supercharges, we conclude that any finite-gap periodic sys- 
tem with even smooth potential is characterized by a hidden bosonized nonlinear N = 2 
supersymmetry of order 2n + 1 [32, 33, 39], 

{Z a , Z h } = 25 ah P 2n+1 (H), a,b = l, 2. (2.16) 

3 Tri-supersymmetric extensions of finite-gap systems 

In this section we show that the application of a non-singular Crum-Darboux transformation 
to a finite-gap system produces a partner system with identical spectrum, and study a 
peculiar supersymmetry appearing in the obtained isospectral pair. 

3.1 Darboux transformations and supersymmetry 

A usual model of supersymmetric quantum mechanics is based on a Darboux transformation, 
by which an (almost) isospectral system can be associated with a given quantum system. 

Consider a Hamiltonian H = — +u(x), and an eigenstate corresponding to a fixed 
eigenvalue E+, Hip* = E+ip*. Here we do not assume any regularity conditions for ■0*. It can 
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be a physical eigenstate, or a second, non-physical solution of the second order differential 
equation, corresponding to a physical energy level E+, or can be a solution corresponding to 
a nonphysical value E*. The Darboux transformation is generated by a first order differential 
operator A 1 = — (In ■?/>*)', which annihilates ^.iV'* = 0, and relates H with another 
Hamiltonian 

d 2 d 2 
H=-— + u(x), u(x) = u(x) - 2— In^, (3.1) 

by means of an intertwining relation 

A X E = HA 1 . (3.2) 

Then for eigenstates of two Hamiltonians corresponding to the same arbitrary value of energy 
E, we have 

HtpE = E^j e , H4> e = E4> E , (3.3) 
$e = ^L— A^e, ij E = -j=^==A\i) E . (3.4) 

Relations (3.3), (3.4) have a symmetry H <-> H , ipE ipE, A <-> A* . This reflects a property 
that the transformation corresponding to the adjoint intertwining relation 

A\H = HA\ (3.5) 

is generated by the operator A\, which annihilates a state i/j* = l/ip*, A\i]j± = 0, and acts 
in an opposite direction by relating H with H . It is easy to see that the both Hamiltonians 
can be represented in terms of operators A\ and A[, H — A\Ai + E*, H = A\A\ + E+. 

Usually, the Darboux transformation is chosen to annihilate a nodeless physical ground 
state ipQ with energy E . In such a case, the potentials u(x) and u(x) are both smooth and 
regular, or both have the same singularities 5 . In a non-periodic case, the physical nodeless 
ground state ipo vanishes at the ends of a (possibly infinite) interval. As a consequence, there 
is no physical partner state with the same energy in the spectrum of H . Indeed, the state 
■00 — l/V'o annihilated by A\ is divergent at infinity and is not physical. In this case both 
systems are almost isospectral, their spectrum is the same except the energy level Eq to be 
absent from the spectrum of H. Note that from the viewpoint of the adjoint intertwining 
relation (3.5), the transformation from H to H is generated by the operator A\ associated 
with a nonphysical state ipo = 1/ipo, which corresponds to a nonphysical for H eigenvalue 
Eq. On the other hand, in correspondence with (3.4), for E = Eq we still have relations 
V>o = A[fj and ^ = l/V'o = A iVo, but fj = f x if)%(x)dx, r] = f 4~ 2 ^ are tne 
non-physical, non-normalizable solutions of the equations Hi] = E r] and Hfj = E fj . 

In the periodic case with ip corresponding to the singlet band-edge state of the lowest 
energy, f] and fj are the non-physical, non-periodic divergent solutions. From this discussion 
it is also clear that if the Darboux transformation is realized with a nodeless state such 
that both states ip± and are not physical (non-normalizable), the energy level E+ is 
absent from the spectra of both partner systems, and physical energy levels satisfy a relation 
E > E±. 



'Singular Darboux transformations generated by the states with nodes also find some applications, see [41]. 
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The relation between the Darboux transformation and the usual supersymmetric quan- 
tum mechanics is direct. The Hamiltonians H and H shifted for the constant E* are known 
as superpartner Hamiltonians, and form a superextended system described by the matrix 
Hamiltonian 

' H+ 

where 



H_ = H-E<, = —^ + W 2 -W', H + = H-E* = —f- + W 2 + W, (3.7) 

ax 2 ax 2 

and W(x) is a superpotential, W = —4- hi if;*. 

With the Darboux transformation, two Hermitian linear differential matrix operators 

Qi = ( Jt A q ) > ^ = i<r 3 Qi (3.8) 

are associated, in terms of which intertwining relations (3.2) and (3.5) take a form of con- 
servation laws for supercharges Q a , [Q a , H] — 0, a — 1, 2. Together with Hamiltonian (3.6) 
they generate the linear N = 2 superalgebra 

[Q a ,H] = 0, {Q a ,Q b } = 25 ab H. (3.9) 

A diagonal Pauli matrix o" 3 plays here the role of the grading operator, [<t 3 , H] = 0, {0-3, Q a } = 
0. 

In a non-periodic case, if if)* or l/if)+ is normalizable, there exists a two-component phys- 
ical state annihilated by both matrix supercharges, which is a ground state of zero energy 
of one of the super-partner subsystems. It is invariant under corresponding supersymmetry 
transformations generated by Q a , and we have the case of exact, unbroken supersymmetry. 
The supersymmetric doublets of states corresponding to positive energies are mutually trans- 
formed by supercharges Q a in correspondence with (3.4). In the case if both if>* and l/if>* are 
not physical, all the states of supersymmetric system (3.6) are organized in supersymmetric 
doublets, including the states of the lowest energy, that takes here a nonzero, positive value. 
This picture corresponds to the broken supersymmetry, which describes a pair of completely 
isospectral super-partner systems. 

In the case of a periodic quantum system with a smooth potential, the supersymmetric 
system (3.6) constructed on the base of the Darboux transformation with a nodeless singlet 
band-edge state will be characterized by zero energy doublet of the states given by the 
columns (0, ^o)' and (1/^0,0)*. Both these states are annihilated by the supercharges Q a , 
and corresponding N = 2 supersymmetry is unbroken. Here the super-partner systems are 
completely isospectral as in the non-periodic case with broken supersymmetry. 

3.2 Higher-order Crum-Darboux transformations 

(Almost) isospectral systems can also be related by differential operators of higher order, 
that corresponds to the situation well described by the generalization of the Darboux trans- 
formation due to Crum. 
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Let a differential operator A k of order k, A k = D k + Y^j=\ cfD k ~ j , annihilates a space 
V spanned by k eigenstates of the Hamiltonian H, V — span{ip 1 , . . . ,ip k }, which are not 
obligatory to be physical. Then, there holds the relation 

A k H = HA k , H = H + 2(cf)' = H- 2(ln W(^, . . . , ^ k ))" . (3.10) 

In the case of the Darboux transformation {k — 1), the Wronskian of a single function is the 
function itself, and (3.10) reduces to the intertwining relation of the standard supersymmetry. 
The operators A k and A\ produce the relations of the form (3.4) for energies E ^ E iy 
i = 1, For k > 1, Eq. (3.10) underlies a higher-order (nonlinear) generalization of 

super symmetric quantum mechanics, see [31, 32, 33]. In a generic case the spectra of H and H 
are almost identical, their spectra can be different in k or less number of physical eigenvalues. 
For a quantum system described by H, one can obtain various partner Hamiltonians H, by 
choosing different sets of eigenstates ipi, . . . , ip k . However, if we want to get the associated 
partner Hamiltonian H with the same regularity properties as the initial Hamiltonian H, 
these states have to be chosen in a special way. 

The higher-order Crum-Darboux transformations can be factorized into the consecutive 
chain of the first order Darboux, or the second-order Crum transformations, see [42] . We are 
interested in the conditions which not only ensure the regularity of the transformations for a 
periodic even finite-gap system, but also produce an isospectral even partner potential. In a 
periodic system, a regular transformation of the second order can be obtained if the kernel of 
the operator A 2 consists of the states corresponding to the edges of the same prohibited band 
[43] . These two states, due to the oscillation theorem mentioned in the previous section, have 
the same period and the same number of alternating nodes. The Wronskian of the functions 
selected in this way is a function of a definite sign, not taking the zero value. Consider a 
Crum-Darboux transformation that annihilates the indicated pairs of the edge-states. In 
the case if the order of the transformation is odd, it has also to annihilate the nodeless 
ground state This guarantees a smooth and singularity-free potential of the partner 
Hamiltonian. 

Concluding, we can construct a whole family of the partner finite-gap periodic systems 
by means of the Crum-Darboux transformations, just by choosing appropriately the singlet 
states of the original system, respecting the rules described above. For instance, the gen- 
erator of a hidden bosonized supersymmetry Z produces a Crum-Darboux transformation 
associated with the trivial selection: it annihilates all the singlets. Since the Wronskian of 
all the singlet states is a nontrivial constant, see Eq. (2.13), we find that the partner Hamil- 
tonian coincides with the original one, and the intertwining relation reduces to the relation 
of commutation of Z and H . 

Below, by means of nontrivial Crum-Darboux transformations, we shall construct a family 
of 2™ — 1 different completely isospectral partner systems for a given arbitrary n-gap periodic 
system, and reveal a special nonlinear supersymmetry appearing in any pair of the total 
family of 2 n systems. The key role in the construction will belong to the already described 
hidden bosonized supersymmetry. 

It is worth to notice here that a regular Crum-Darboux transformation for finite-gap 
periodic systems can also be produced by making use of certain Bloch functions [28, 29, 
44, 45, 46]. A partner for parity-even potential obtained in such a way in a generic case, 
however, will not be an even function. We shall return to this point in the last section. 
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3.3 Tri-supersymmetric extensions and centrally extended nonlin- 
ear N = 4 super symmetry 

Consider an n-gap periodic system, and mark r < n prohibited bands in its spectrum. 2r 
singlet physical states at the edges of these prohibited bands span a 2r-dimensional linear 
vector space which we denote by V + . The Wronskian of the corresponding 2r singlet band- 
edge states is a nodeless 2L-periodic even function. Let Q + be a linear differential operator 
of order 2r that annihilates the space V+, 

2r 

Q + = D 2r + Y,cf(x)D 2r ->, Q + V + = 0. (3.11) 

3=1 

Singlet band-edge states are periodic or anti-periodic, and can be presented by real wave 
functions. So, the coefficient functions in (3.11) are real. Taking also into account that any 
band-edge state has a definite parity, one can show that (3.11) is a 2L-periodic even operator, 
[T, Q + ] = [R, Q + ] = 0, see Appendix A. 

The kernel of the integral Z has a form Ker Z = V + © V-, where V_ is a supplementary 
2(n — r) + 1-dimensional linear vector space spanned by the rest of singlet band-edge states. 
Then Z can be decomposed as Z = S*Q + , where S' is a differential operator of order 
2(n-r) + l with the property S*Q+V_ = 0. Hermiticity of Z and Eqs. (2.10), (2.15), (3.11) 
mean that 

Z = S*Q + = Q{S, (3.12) 

and 

2(ra-r)+l 

-iS = D 2{n - r)+1 + c s j {x)D 2{ - n -^ +1 -\ (3.13) 

3=1 

where the coefficient functions are real and cf(x) = cf(x). From the properties of Z and 
Q + , we also find that S is a 2L-periodic parity-odd operator. 

Now we show that KerS* = V_. To this end we note that in accordance with Eq. (3.10) 
and the Abel identity W = —ci(x)W, the equality cf(x) = c^{x) obtained directly from 
(3.12) means that the application of the Crum-Darboux transformations with the operators 
Q + and S produces the same non-singular super-partner Hamiltonian H = H + 2(c^)' 
satisfying the intertwining relations 

Q+H = HQ + , SH = HS, (3.14) 

Q{H = HQ\, S ] H = HS*. (3.15) 

The Hamiltonian H describes a periodic system with an even potential of the period 2L 
with n gaps in the spectrum. Thus, there exists an odd Hermitian differential operator Z 
of the form (2.10) commuting with H. The intertwining relations (3.14), (3.15) provide an 
alternative, two-term decomposed form for Z. Indeed, we get [H, SQ + ] = [H, Q+S*] = 0. 
Both operators SQ+ and Q + S^ are of order 2n + 1, and should coincide with Z up to some 
polynomial in H . However, they anticommute with reflection operator R, and this implies 
that 

Z = SQl = Q+SK (3.16) 



11 



Take any of 2(n — r) + 1 singlet states G V- not annihilated by Q + . Multiplying the 
equation = by from the left and using (3.12) and (3.16), we get SQ+Q + ^5/^ = 0. 

The operator Q+Q+ is 2L-periodic operator, and as it follows from (3.14), (3.15), commutes 
with the Hamiltonian H. It changes neither energy nor the period of a singlet state and 
then Q\Q + ^>- = where a is a non-zero number. Hence, SQ+Q + ^- = aS^fr- = 0, 

and we conclude that Ker S = V-. Changing the notation, Q_ = S, we have 

Z = Q\Q- = QlQ+, Z = Q + Ql = Q_Q{, (3.17) 

and 

Ker Q + © Ker Q_ = Ker Z, Ker Q\ © Ker Ql = Ker Z. (3.18) 

This result means a complete isospectrality of the finite-gap periodic systems described by 
the Hamiltonians H and H . Indeed, in accordance with the properties of the Crum-Darboux 
transformations, the action of both operators Q + and Q_ on any doublet of eigenstates of H 
from the interior of permitted bands transforms it into a doublet of eigenstates of H with the 
same energy value. The adjoint operators Q+ an d Q- act in the same way in the opposite 
direction. The singlet states of H annihilated by Q + (or QJ) are transformed by Q- (or 
Q + ) into zero modes of Q\ (or Q\_) being the singlet states of H of the same energy. The 
same picture is valid for the singlet states of H annihilated by Q\ (or Q_) and transformed 
by Q_ (or Q + ) into the corresponding singlet states of H. 

The intertwining relations (3.14), (3.15) as well as factorization of the supercharges of 
bosonized supersymmetry can be rewritten in a compact form once we use the matrix for- 
malism and define an extended Hamiltonain TC and operators Q± and Z, 

«-(!£)■ *-(',!)■ <M« 

Here relations (3.14), (3.15) and (3.17) can be presented as 

[H,Z] = 0, [H,Q±} = 0, (3.20) 

Z=Q_Q+ = Q+Q- (3.21) 

The triplet Q + , Q_ and Z is a set of commuting integrals for the superextended system 
described by the matrix Hamiltonian 7i. There exists a common basis, in which Q±, Z and 
H are diagonal, and since all these operators are self-adjoint, their eigenvalues are real. We 
have a chain of equalities 

2n+l 

Z 2 = Q+Q-Q+Q- = Q 2 + Q 2 _ = P Z {H) = JJ (H - Ej), (3.22) 

3=1 

where Pz is a positive-semidefinite spectral polynomial, and Ej are energies of the singlet 
states of subsystems. In correspondence with (3.22), the band-edge states of the extended 
system are organized in supersymmetric doublets, on which Z takes zero values. The states 
of the interior of permitted bands are organized in energy quadruplets, on certain pairs 
of which Z takes nonzero values ±y / P z (E). The diagonal components of Q 2 ± consist of 
differential operators of orders 4r and 4(n — r) + 2. One of these two numbers is less than 
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2n + 1. Suppose that it is the case of the operator Q\. Its lower diagonal component 
Q\Q + satisfies a relation [H, Q\Q + ] = 0. According to the general theory of finite-gap 
systems, the only operators commuting with H of order lower then 2n + 1 are polynomials 
in the Hamiltonian, and we conclude that Q + Q + is such a polynomial, which has to take 
zero values on 2r singlets belonging to KerQ+. The same arguments hold for the upper 
component of the operator Q 2 +) and we find that 

2r 

Ql = P + (7i) = \[(H-E+). (3.23) 
i=i 

With making use of (3.22) and (3.23), we get also 

2(ra-r)+l 

Ql = P_(r L )= J! (H-E-), (3.24) 

3=1 

where P±{7i) are positive-semidefinite operators, and E ± are the energies of the correspond- 
ing band-edge states annihilated by Q±. The eigenvalues of Q± are ±^P±(E), where the 
signs of square roots are correlated with the square root sign of corresponding eigenvalue of 
Z in accordance with Eq. (3.21). 

In addition to the non-trivial integrals of motion Z and Q±, the Hamiltonian H possesses 
another triplet of trivial, mutually commuting, integrals 

ri = <7 3 , T 2 = R, T 3 = Ra 3 , (3.25) 

which satisfy the relations Tf = 1, i = 1,2,3, Tilers = 1. Any of r« can be chosen as 
a Z2-grading operator r*. Any of the non-trivial integrals of motion either commutes or 
anti-commutes with any of the trivial integrals. Fixing the grading operator, we classify any 
nontrivial integral as a bosonic or fermionic, while the Hamiltonian and the trivial integrals 
(3.25) are always identified as bosonic operators. In correspondence with this identification, 
a certain superalgebra is generated. For all the three possible choices of the grading operator, 
one of the nontrivial integrals plays the role of a bosonic, Z 2 -even operator, while the other 
two integrals are classified as fermionic, Z 2 -odd operators, see the Table below. We name 
this structure the tri- super symmetry 6 . 



Grading operator 


<?3 


R 


a 3 R 


Bosonic integral 


Z 




Q- 


Fermionic integrals 


Q+,Q- 


Z, Q- 


2, Q+ 



The complete structure of the tri-supersymmetry will be described in the next subsection. 
Here, let us choose = cr 3 as the grading operator, and discuss the corresponding nonlinear 

6 Such a structure was observed for the first time in the N — 2 superextended Dirac delta potential 
problem [47], where the basic triplet of nontrivial integrals has a completely different nature. In particular, 
there both supercharges of the hidden bosonized N = 2 linear supersymmetry of the form (3.9) are nonlocal 
operators, cf. (2.16) and (2.15) 



13 



super symmetric subalgebra generated by the local integrals of motion, forgetting for the 
moment the nonlocal integral R. Introduce the notation 

Q { ± = Q±, Ql ] =i^Q ± . (3.26) 

These fermionic supercharges together with bosonic operators Z and TL generate the super- 
algebra 

{Q+\ Q { 1 ] } = 25 ab P + (H), {Q { *\ QL 6) } = 25 ab P4H), (3.27) 
{Q+\ Q { - ] ] = 25 ab Z, (3.28) 
[H, Q ( ± ] ] = [H, Z] = [Z, Qi a) ] = 0. (3.29) 

Superalgebra (3.27), (3.28), (3.29) is identified as a centrally extended nonlinear N = 4 
supersymmetry, in which Z plays a role of the bosonic central charge 7 . 

The supercharges annihilate a part of the band-edge states organized in supersym- 
metric doublets, while another part of supersymmetric doublets is annihilated by the super- 
charges q!l\ The band-edge states which do not belong to the kernel of the supercharges 
<2+^ (or <2^) are transformed (rotated) by these supercharges within the corresponding su- 
persymmetric doublet. The bosonic central charge Z annihilates all the band-edge states. 
So, we have here the picture reminiscent somehow of the partial supersymmetry breaking 
appearing in supersymmetric field theories with BPS-monopoles [49]. 



3.4 Tri-supersymmetry and su (2\2) 

Let us study the complete algebraic structure of the tri-supersymmetry. To do this, consider 
the set of integrals of motion created by the multiplicative combinations of the trivial integrals 
with nontrivial ones, 

H, Vi, T a Z, T a Q+, (3.30) 

where a = 0, 1, 2, 3, and by T we denote a unit two-dimensional matrix. Each of these 
integrals either commutes or anticommutes with any of Tj defined in (3.25). Identifying one 
of Tj as the Z 2 -grading operator r*, we separate the set (3.30) into eight Z 2 -even (bosonic) 
operators commuting with and eight Z2-odd (fermionic) operators, which anticommute 
with r*. Though this separation depends on the choice of the superalgebra in all three 
cases has, in fact, the same structure. 

To reveal this common superalgebraic structure for all the three possible choices of the 
grading operator, we denote the fermionic operators as (Fx,..., F 8 ). The set of bosonic 
operators we write as (7i, r*, Si, £2, Bi, ... , B4), where £^2 are two trivial integrals from 
the set (3.25) to be different from r*. Denote by Pb(J~C) a universal polynomial produced 
by the square of any of the four integrals B a , B\ = Pb(H), a = 1, ... ,4. Given T*, we 
can separate the set of fermionic operators in two subsets depending on the commutation 

7 The basic structure of the algebra of finite order differential operators of a form more general than 
(3.27)-(3.29) was discussed by Andrianov and Sokolov [48], but outside the context of finite-gap periodic 
systems and parity-even potentials. In comparison with (3.27)-(3.29), the superalgebraic structure of Ref. 
[48] includes some additional independent polynomial of the Hamiltonian. As a consequence, instead of the 
relation Z 2 = P + (H)P-(H), which follows from Eqs. (3.22)-(3.24) and reflects the nature and peculiarities 
of the band structure [34], its analog in [48] has a different form, see Eqs. (40), (41), (43) and (46) there. 
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relations with the integrals Si and E 2 . The subset which commutes with E 2 , we label as 
Fp, /i — 1, . . . , 4, [E 2 , F^] = 0, and denote a universal polynomial corresponding to a square 
of any of these fermionic operators by P 22 , F£ = P 22 . For the subset of fermionic operators 
which commute with E 1; we put the index F\, A = 5, ... ,8, [Ei,.F\] = 0, and denote the 
analogous universal polynomial by Pn, F% = Pn. Identifying in the described way the 
integrals F M , F\, B a , Ei and E 2 , and computing directly all the (anti) commutators, we find 
the superalgebra of tri-supersymmetry, which may be presented in the same form modulo 
some special polynomials in dependence on the choice of T*. These additional polynomials we 
denote by Pi 2 , P\b and P 2 b, where the subindex with two entries indicates the origin of the 
(anti) commutation relation. Polynomial Pi 2 comes from the anticommutators of F\ with F M , 
polynomial P 2 b comes from the commutators between the integrals F M , ([E 2 ,F M ] = 0) and 
B a , polynomial P\b has analogous sense. Since the polynomials Pi 2 , P\b and P 2 b depend on 
the grading, their explicit form together with explicit form of bosonic and fermionic operators 
for all three choices of the grading operator are presented in Appendix B, see Tables 1, 2 and 
3. With the described notations, the anti-commutation relations for fermionic operators, 
and commutation relations between bosonic and fermionic operators are presented in Tables 
4 and 5. 

The identification of the complete superalgebra of the tri-supersymmetry can be achieved 
now if we analyze the still missing commutation relations between Z 2 -even generators. In- 
troduce the following linear combinations of them, 

Q? ] = \(Bi± B 3 ) = \BJl±, = ~ (B 2 ± P 4 ) = ~B 2 n±, (3.31) 

4 ±] = \ (Si ± S 2 ) = isill ± , (3.32) 
where H± = |(1 ± r*) are the projectors. These operators satisfy the following algebra 

= i4 ±) P B (H), (3.33) 
= ie bGl ± \ a, b =1,2, (3.34) 

[G ( a + \ Gt } ] = [4 + \SP] = [Jt\ G ( a +) ] = 4^] = 0, (3.35) 

The commutation relations (3.33)-(3.35) correspond to the direct sum of two deformed su(2) 
algebras in which 7i plays a role of a multiplicative central charge. This bosonic subalgebra 
is reminiscent of the nonlinear algebra satisfied by the Laplace-Runge-Lenz and angular 
momentum vectors in the quantum Kepler problem [50, 51]. 

It is known that in the case of the quantum Kepler problem, its nonlinear symmetry 
algebra is reduced on the subspaces of fixed energy E < 0, E = and E > to the Lie 
algebras so(4), so(3, 1) and e(3), respectively, where e(3) is the 3D Euclidean algebra. Let us 
see what happens with our tri-supersymmetry under similar reduction. First, consider any 
4-fold degenerate energy level E ^ Ei corresponding to the interior part of any permitted 
band. Rescaling the operators, — > = /P B (E), we find that together with 

they generate the Lie algebra su(2) © su(2). These operators satisfy the relations 
Ji + '' ] Ji +S> — fn+, J-~^J-~^ = |n_, where the summation in % — 1,2,3 is assumed. The two 



7 (±) c(±) 
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common eigenstates of the Hamiltonian, with energy E ^ Ei, and of the grading operator 
r*, with eigenvalue +1 or —1, carry the 1/2 © 0, or © 1/2 representations of su(2) © su(2), 
where the first (second) term corresponds to the generators (J^). The fermionic 

generators mutually transform the states from the two eigenspaces of the grading operator. In 
accordance with the total number of independent fermionic generators, the energy subspace 
with E ^ Ei carries an irreducible representation of the su(2|2) superunitary symmetry, 
which is a supersymmetric extension of the bosonic symmetry u(l) © su(2) © su(2), where 
the u(l) subalgebra is generated by the grading operator, see Ref. [52]. Having in mind that 
the Hamiltonian appears in a generic form of the superalgebra as a multiplicative central 
charge, we conclude that the system possesses a nonlinear su(2|2) superunitary symmetry 
in the sense of Refs. [31, 32, 51]. 

If we reduce our extended system to the subspace corresponding to any doubly degener- 
ate energy level Ei corresponding to a doublet of band-edge states, the bosonic part of the 
superalgebra is reduced to the algebra u(l) © e(2) © e(2), where the first term corresponds 
to the integral r*, while other two correspond to the two copies of the 2D Euclidean alge- 
bras generated in the eigensubspaces of by the rotation operators and commuting 
translation generators ■ Note that the supersymmetry of the form similar to the present 
one reduced to a level Ei was analyzed in [53] in the context of spontaneous supersymmetry 
breaking in 3+1 dimensions. 

3.5 Self-isospectrality conjecture 

In the realm of supersymmetric quantum mechanics associated with linear superalgebraic 
structure, the complete isospectrality in the non-periodic systems is related to the super- 
symmetry breaking, that means that the doublet of the ground states is not annihilated by 
supercharges. In [25] , Dunne and Feinberg considered supersymmetric extensions of periodic 
potentials. They argued that in contrary to the usual situation, the complete isospectrality of 
super-partner Hamiltonians could appear without violation of the supersymmetry. As one of 
the examples of the situation they presented a one-gap Lame Hamiltonian, where the super- 
symmetric extension was provided by the first-order Darboux transformation corresponding 
to the first order supercharges defined by (3.26). The super-partner Hamiltonian showed 
to be the original one but displaced for a half of the period. As we mentioned at the be- 
ginning, such a phenomenon of a half-period displacement of super-partners was named in 
[25] the self-isospectrality. We explained above how the complete isospectrality emerges 
due to the tri-supersymmetry, namely, its local part (3.27)-(3.29). In this framework, the 
symmetries of the N = 2 superextended one-gap Lame system have to be completed by 
adding two other supercharges of order 2 and a bosonic integral Z of order 3. The 
second order supercharges do not annihilate the doublet of the ground states, and the 
tri-supersymmetry is spontaneously partially broken. As we showed, this turns out to be 
a general feature of the tri-supersymmetric systems constucted by extension of a finite-gap 
periodic system by means of a regular Crum-Darboux transformation. 

We could ask, motivated by [25], for the indications on the self-isospectrality in the 
tri-supersymmetric extensions of the finite-gap systems. In our current setting, the self- 
isospectrality arises if the translation in the half-period L provokes inversion of the Wron- 
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skian, 

WH* + L)=C± w ^ y (3.36) 

where C ± are some nonzero constants. Indeed, such a displacement produces changes in 
the sign for coefficient functions cf(x) = — (lniy ± ) / of the operators Q±, and therefore 
transforms the latters into their conjugates, 

cf(x + L) = -cf(x), Q±(x + L) =Ql(x). (3.37) 

Making the translation for L in the intertwining relations and comparing the result with 
their conjugates, we reveal that 

H(x) = H(x + L), (3.38) 

and therefore the self-isospectrality does appear. 

The construction considered in this section provides a receipt how to get isospectral tri- 
super-symmetric partners for a given n-gap Hamiltonian. We have seen how the partner 
Hamiltonian H is determined uniquely once we make a separation of the singlet states into 
two disjoint families. There exist Ylk=o (fe) = 2 n distinct separations which respect the 
rules explained in the subsection on the Darboux-Crum transformations. Since one of them 
is trivial (includes all the singlet states and corresponding integral Z commutes with the 
Hamiltonian H) we end up with 2 n — 1 tri-supersymmetric isospectral extensions of the 
given n-gap system. All the isospectral extensions can be obtained by successive first order 
Darboux and second-order Crum transformations. 

If antiperiodic singlet states are present in the spectrum, among the possible separations 
of the singlet states there exists an exceptional one, given by sorting out the singlets into 
mutually orthogonal families of periodic and antiperiodic states. Despite the lack of the 
proof, we conjecture that this "natural" separation leads to the self-iso spectral supersymmetry 
characterized by the partner Hamiltonian H to be the original one but displaced in the half 
of the period. 

Suppose an n-gap system H with n > 1 has antiperiodic singlet states in its spectrum, 
and H is a shifted for the half-period Hamiltonian obtainable by the Darboux-Crum trans- 
formation associated with the specified natural separation of the singlets. Let Q± be the 
generators of the Crum-Darboux transformation associated with a separation of singlets dif- 
ferent from the natural one. Then shifted for the half-period operators Q± will generate a 
corresponding Crum-Darboux transformation for the system H . In such a way we obtain a 
new pair of self-isospectral systems Hq and Hq: 



x^x+L 




Including H and H, we can get 2 n_1 distinct pairs of self-isospectral Hamiltonians. Starting 
with any n-gap system system that has nonzero number of antiperiodic singlet states, we 
would be able to construct 2 n_1 extended self-isospectral Hamitonians H. 
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In the next section, our self-isospectrality conjecture will be supported by the study of 
the tri-supersymmetry of the associated Lame equation. 



4 Associated Lame equation and its isospectral exten- 
siones 

We apply here a general theory developed in the previous sections to a broad class of finite-gap 
systems described by the associated Lame equation. In particular we study the isospectral 
extension based on the natural separation of the singlet states into periodic and anti-periodic 
ones, and show that it leads to the self-isospectral tri-supersymmetric systems. We provide 
an explicit form of both the diagonal and non-diagonal integrals of motion of the extended 
system. The examples of isospectral extensions not possessing a property of self-isospectrality 
are presented as well. 

Associated Lame equation is a two-parametric second order differential equation of Fuch- 
sian type with four singularities and doubly-periodic coefficients, 



where C m = m(m + 1), C\ = 1(1 + 1) are real numbers and dnx = dn (x, k) is Jacobi 
elliptic function with modular parameter k G (0, 1); k' G (0, 1) is a complementary modular 
parameter, k' 2 = 1 — k 2 . Lame equation (1 = case), obtained originally by separation of 
the Laplace equation in elliptical coordinates, has been a subject of extensive studies with 
use of both analytical [54, 55] and algebraical [8] methods. Due to appealing properties 
of its solutions, the equation (4.1) found the applications in diverse areas of physics. In 
solid-state physics [7], it represents a stationary Schrodinger equation of a model of one- 
dimensional crystal with a more realistic potential than Kronig-Penney or Scarf potentials. 
This equation, especially its I — case, plays an important role in many other fields of 
physics as well. For instance, it appeared in some expansions of scattering amplitudes [56], 
in the study of bifurcations in chaotic hamiltonian systems [17], it governs distance red-shift 
for partially filled-beam optics in pressure-free FLRW cosmology [57], it was used in the 
study of static SU (2) BPS monopoles [9] and kink solutions [58] in the field theory. 

4.1 Construction of self-isospectral extension 

The spectrum of the one-dimensional periodic system governed by the Hamiltonian operator 
corresponding to (4.1) 



consists of the valence bands and the prohibited zones (gaps) which alternate mutually until 
energy reaches a semi-infinite band of conductance. Configuration of the spectral bands 
depends sensitively on the constant parameters. As long as m and I acquire integer values, 
which we suppose to be the case from now on, the spectral bands are arranged such that 
only finite number of gaps appear. The period 2L of the potential with C m ^ Ci in (4.2) 
is equal to 2K, where K = f^ 2 (l — k 2 sin 2 (p)~^d(f) is the complete elliptic integral of the 




(4.1) 



H m ,i = -D 2 - C m dn 2 x - Q 



k' 2 



(4.2) 



dn 2 x 
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first kind. The case C m = C\ corresponds to the Lame system with the same value of C m 
but C\ — and the period 2L = K; it is discussed separately in the Appendix C. The 
independent change of parameters m — > — m — 1, I — > — I — 1 leaves the Hamiltonian (4.2) 
invariant so that we can consider the case m > I > without the loss of generality. In this 
case the system is m-gap. 

To start on the construction of the tri-supersymmetric extension, we focus to the sep- 
arations of the band-edge states. As we announced, the separation into the periodic and 
anti-periodic singlets will be used here. Construction of the operators Q + and Q- associated 
with any separation would require an explicit knowledge of the band-edge states, i.e. an 
explicit solution of the stationary Schrodinger equation. Fortunately, in the case of natural 
separation this rather compelling work can be passed with the use of peculiar properties of 
the model. 

The present one-dimensional system is closely related to the finite-dimensional represen- 
tations of Lie algebra sZ(2,R). The Hamiltonian (4.2) can be written as a second order 
polynomial in generators of a finite-dimensional irreducible representation of sl(2,M). This 
important feature underlies quasi-exact solvability of the model, implying that a finite num- 
ber of eigenstates corresponding to band edges can be found by purely algebraic means 
[59, 60]. For integer values of m and I, m > I > 0, the space of 2m + 1 singlet states of 
the associated Lame system can be treated as a direct sum of two irreducible non-unitary 
representations of sl(2, M) algebra of dimensions m — I (spin j_ = \{m — I — 1)) and m + l + 1 
(spin j + = \{m + l)) [61,62]. 

This fact is deeply related to the structure of the band-edge wave functions of the system, 
m + I + 1 of which can be factorized formally as 

*„ = ^(0, ( 4 -3) 
whereas the remaining m — I singlets acquire the form 

= ^(0- (4-4) 
Here we introduced the functions a = ""I' 1 , v = cn m_/_1 xdn l+1 x and a new variable 

' dm 1 

£(x) = that varies smoothly from — oo to +oo in the period interval (—K,K). The 
functions an d ^v(0 are, in general, polynomials of order m + I and m — I — 1 in £ 

and, as we will see, lie in vector spaces of the irreducible representations of s/(2,M) of the 
dimensions m + l + 1 and m — I. 

As £ = £(x) is periodic, the factors fi and v dictate periodicity or antiperiodicity of the 
eigenfunctions; for even m + I the wavefunctions (4.3) are periodic while the functions (4.4) 
are anti-periodic. For m > I > 0, function fi has one node in the interval (— K, K), while 
v can have at most one node there. Being the polynomial in £, the function (T v ) can 
acquire at most m + l (m — l — 1) zeros in this interval. Combine these facts with the general 
properties of periodic and anti-periodic states of Hill's equation discussed in Section 2. The 
resulting picture shows that starting with periodic ground state and anti-periodic states at 
the edges of the first gap, the gaps with periodic and anti-periodic states at their edges alter 
with energy increasing until the (m — /)th gap is reached. For higher gaps, all the remaining 
singlets at the edges are of the same nature as edge states of the (m — /)th gap, periodic or 
anti-periodic with even or odd number of nodes, see also ref. [34]. 
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To reveal the algebraic form of (4.2), we shall recover how the Hamiltonian acts on the 
"dynamical" part (JF V ) of the wavef unctions. Transforming out the function \i and writing 
the result in the variable £ , we obtain 

h, = ^)- 1 H m ^ = -{k>\T+) 2 + (l + k' 2 )(T ) 2 + (T-) 2 + k\l-m)T ) 

= - (1 + kX) (1 + £ 2 ) ^ + m' 2 (m + l - l)e + 2(to + Ik' 2 - 1) + A; 2 )| 

-A;' 2 (m + /)(m + /- 1)£ 2 + const. (4.5) 
Here the operators 

T+ = ^-(m + Z)£, T° = ^-^±i, T- = ^, (4.6) 

[T + ,T~] = -2T°, [T°,T ± ]=±T ± , (4.7) 

are the generators of irreducible representation of sl(2, M) acting on the vector space spanned 
by monoms {1,£, ..,£ m+ '}. Representation (4.6) is specified by the eigenvalue + 1) of 

the Casimir C = — (T ) 2 + \(T + T- + T~T+) corresponding to sZ(2,R) spin j + = \{m + /). 

The other algebraic form of (4.2), h v , that acts on the "dynamical" part T v of the wave 
functions (4.4), can be obtained by performing the gauge transformation with the other 
common factor v. Alternatively, we can use the apparent symmetry /x|;_ > _;_ 1 = v and write 
down h u immediately just by substituting I — > — I — 1 in (4.5) and (4.7), 

K = h,\ M ^ = {u)- l H- hl u = - (k' 2 (T + ) 2 + (1 + A;' 2 )(T ) 2 + (f-) 2 + k 2 (-l - 1 - m)T°) . 

(4.8) 

We denoted by T p = T p |^_2_i, p = 0, +, — , the s/(2,R) generators of (m — I) -dimensional 
representation, where T p are the generators (4.6) of spin-j + representation. Note that the 
"effective" Hamiltonian is Hermitian with respect to a scalar product defined with a 
nontrivial weight, (f,g) = /* (OsCfX 1 + k' 2 f)- l+ ^{\ + £ 2 )~ m +^£ ; the same is true for 
h u with the change I — > — I — 1. 

Now the background of the natural separation of the singlet states is clear. The periodic 
and antiperiodic singlet states carry two different irreducible representations of s/(2,R) of 
dimensions m + 1 + 1 and m — l. The number of the periodic and anti-periodic singlet states 
depends on the values of m and /, while their total number 2m + 1 is fixed by the number 
of gaps to. For instance, for to = 3, 1 = there are to — I = 3 periodic band-edge states 
while for to = 3, / = 1 we have to + I + 1 =5 periodic singlet states. To avoid possible 
confusions during the construction of the supercharges, let us change the notation slightly. 
We will denote by X~ t an operator which annihilates all the functions (4.4), and the operator 
annihilating all the states (4.3) will be Y~ v 

First, let us consider eigenstates (4.3) covered by the m + l + 1 dimensional representation 
of s/(2,R). An operator of the order m + l + 1 which annihilates the representation space 
spanned by the monoms {1, £, . . . , £ m+z } has the following general form 

y~ m ,i = <*m,id? +l+1 . (4.9) 

The function a my i is fixed uniquely as we require the coefficient at ]j m+l+1 of the operator 
Y m,i = i"2/m,(J |f=c(x) to be equal to one. It reads explicitly a myi = (^f) m+ ' +1 - We present 
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below two equivalent forms of the operator Y~ v The second, factorized expression, will be 
particularly helpful in study of the limit case k — > 1. An explication how to get it from (4.9) 
can be found in [63], 



Hn m+1 r /mV \ m+l+1 rWr 
v - _ n m+i+i , y n m+i+i-j _ an__x_ / cn_x \ (in.) 



cn m+i x 



(m+0/2 



n/A f k 2 (m — l)cn 2 x ,,, 9 , 9 A sna; \ 
V° + (^V J(i +dni) j^)' (4 ' 10) 

The upper index of the ordered product corresponds to the first term on the left side while 
the lower index denotes the last term on the right side of the product. We can construct 
the operator X~ t in the same way or just by making the substitution I — > — I — 1 in (4.10) 
which interchanges considered algebraic schemes. Explicitly, we get 

m-l , m +l / 2 \ m ~ l j -l-l 

an x f cn x \ dn x 



y- - nm-l | \^ „X nm-l-j _ dn? " + x ( cv?x n Y 

m ' 1 ~ + ^ ° 3 ~ cn^'^x V dna; 

j=i x - 



cn m 1 1 x 



n/A f k 2 (m + I + l)cn 2 x /2 2 A snrr \ 

— H ><* +inx) )^^)- (4n » 

J= _( m _i_l)/2 ^ v / / 

As we explained in the preceding section, the coefficients of the second highest derivative of 
X~ l and Y ml coincide and enter the explicit construction of the superpartner Hamiltonian, 

x y W m,i , 2 ( C m ~ Ci) snxcnx 1 2 cnxsnx 
Cl " Cl =Cl ="WW =A; ^ dn^ = 2 (m - /)(m + / + 1)A; Ad^T- (4 ' 12) 

The equality cf = cf reflects the coincidence of the Wronskians of the kernels of X~ ; and 
Y ml up to inessential numerical factor related to the arbitrariness in normalization of their 
zero modes. The essential part of these Wronskians is given by the nodeless function 

W mJ i{x) = (dnx)^ m - l ^ m+l+1 \ (4.13) 

whose invariance with respect to the change I — > — I — 1 just reflects the indicated equality 
of the coefficients. 

The Jacobi function property dn (x+K) = k' /dnx shows that the Wronskian W m> i(x) sat- 
isfies the relation (3.36), and supports our conjecture that the natural separation of the singlet 
states into periodic and anti-periodic states results in the self-isospectral tri-supersymmetric 
system with the partner Hamiltonian operator H = H^ nl to be the original Hamiltonian 
translated for the half-period, H^ t (x) = F mi (i + K). Its explicit form is 

K,i = H~ tl + 2c[ = —D 2 - Qdn 2 x - C m -^-. (4.14) 

dn x 

Recalling (3.17), the generator of the hidden bosonized supersymmetry of the associated 
Lame system (4.2) acquires the following factorized form 
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2_ dn 1 x ( cn 2 x d \ m 1 ( dnx\ 2m+1 / cn 2 x d \ m+l+1 dn l x 

m ' 1 cn m-i+i x \ ^ nx g x J y cn x / V dnx dx I cn m+ ' x 



dn / sn 2 :r d \ m 1 f dnx\ 2m+1 ( sn 2 x d \ m+l+1 d n l x 
sn m-/+i x y ^ y ysnx/ y dnx / sn m+ ' a; ' 



m+l / 2 \ — 2 — 1 

where we used alternative expressions X m i (a;) = s ^Z l+ ? x (j^D) ^-i-i x and Y~^(x) 

m+l ( 2 \ m +'+l ; 

( ^jf D ) sn^+'x ' obtained with the use of a specific identity 



/sn 2 a: V 1 1 /cn 2 x \ J 1 



sn-? +1 a; y dnx / $a?~ x x cn.i +1 x y dnrr ) va?~ x x 

With making use of the same identity we can prove that Y~ t (x + K) = (— l) m+l+1 {Y^ ^x))* 
and X^ ^x + K) = (— l) m ~ l (X m j(x)y . Finally, we write down the obtained extended Hamil- 
tonian 7i as well as the Hermitian diagonal and antidiagonal integrals 

Zm,l = i^( Z ^ ° )=fm + l( Z^X + K) \ (4 _ ig) 



o z- }l j v Z m ,i( x ) 

X* = im ~ l ( v+ X ^)=^( ir-fL ^ X ^ (X) )• (4.20) 



which represent an explicit realization of (3.19). In accordance with the analysis of Section 
3, the integrals Q± are identified with X m l and y m j in the following way: <2_ = X m X and 
Q+ = y m ,i when m — I is odd, and their roles are interchanged when m — I is even. 

4.2 Self-isospectral pairs and "superpotential" 

Let us indicate on an interesting representation of the self-isospectral pairs of the Hamiltoni- 
ans that generalizes a representation H± = —D 2 +W 2 ±W of the super-partner Hamiltonians 
in terms of the superpotential in the case of the usual (linear) N = 2 supersymmetry. 
The self-isospectral pair (4.2), (4.14) can be presented in the equivalent form 

= -° 2 + 2 j^W {lnWm " y2 ± {lnWm ' l) " + + k ' 2) \ {Cm + Clh (4,22) 

where W m> i is the Wronskian (4.13) corresponding to the kernels of operators X ml and Y~ t . 



Let us denote C + = \J\{C m + Cj), C_ = \{m — l)(m + / + 1) = \{C m — C{), and define a 
function 

W = - (\ndn c+ x)' . (4.23) 
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Then (4.22) can be rewritten equivalently as 

H m i ~ (1 + k ' 2 ) C l = -D 2 + W 2 ± ^-W. (4.24) 

Eq. (4.24) is reminiscent of supersymmetric quantum mechanics representation (3.7). This 
is not just a coincidence. In the case m — I = 1, the system H~ _ ± is characterized by 
the presence of only one periodic singlet = dn m x, which is the ground state with energy 
corresponding to a subtracted constant term on the left hand side of Eq. (4.24). In this 
case the first order supercharge X m , m -i reduces to one of the first order supercharges (3.8) 
of N = 2 supersymmetry, (4.23) takes a form of a usual representation of a superpotential 
in terms of the ground state, and (4.24) reduces to (3.7). 

There exists a simple generalization of a classical model for supersymmetric quantum 
mechanics to the case of nonlinear supersymmetry of order n > 1 [32]. It consists in the 
change of the boson-fermion coupling term 9 + 9~W in classical Hamiltonian for nO + 0~W, 
where 9 + 9~ is a classical analog for a 3 and 9^ are Grassmann variables describing fermion 
degrees of freedom. However, unlike the linear case n — 1, for n > 1 such a generalization 
suffers a problem of the quantum anomaly, which can be solved in a general form only for 
n = 2 [33]. One can show that for m — I — 2, when the supercharge X m , m -2 is the differential 
operator of the second order, representation (4.24) is in correspondence with the solution of 
the quantum anomaly problem for n — 2. 

Let us stress that for m — I > 1 the argument of logarithm in the definition of the 
superpotential-like function (4.23) does not correspond to a ground state of the system 8 . 
It is just the appropriately rescaled logarithmic derivative of the Wronskian (4.13), W = 

-g^lnWW)'- 

4.3 Some examples 

Due to the general result presented in Section 3, there exist 2 m — 1 isospectral extensions of 
an m-gap associated Lame Hamiltonian. The new systems can be obtained from the initial 
system by sequent use of Darboux transformation and the transformation of Crum of the 
second order. Moreover, the complete set of 2 m isospectral systems can be sorted out into 
2 m_1 self-isospectral tri-supersymmetric pairs. We present here an example of Lame system 
to illustrate this picture explicitly. 

First, we explain briefly some subtleties related to the sequent use of the transformation 
of Crum. Let us consider a two-gap associated Lame system represented by Hamiltonian 
H . There are three admissible separations of the singlet states: sorting out the states at 
the edges of the first gap, or of the second gap, or of both of them. Let us denote the 
supercharges associated with these separations as Q±j, j = 1,2,3. Then we can construct 
three isospectral super-partner Hamiltonians satisfying 

H U )Q±,j = Q±,jH. (4.25) 

We can repeat the procedure with anyone of the new systems, obtaining another isospectral 
Hamiltonian H^j) (j refers to the system H^, k denotes the next choice of separation). 

8 For the explicit form of the ground states of the associated Lame system with some values of m and I 
see [62]. 
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Although we can make this procedure repeatedly, it will generate only limited number of 
new systems. This is due to simple rules following the sequent use of transformations of 
Crum. 

The first rule: when we choose sequently two identical separations, we return to the initial 
system. In our two-gap setting, let us start with the first separation and construct the new 
Hamiltonian Hi with the help of operator <5+,i- The operator which annihilates the states 
at the edges of the first gap of Hi coincides with Q+ v So, repeating the procedure with Hi 
with the same separation, we obtain a new Hamiltonian Hnu which is related with initial 
H by the intertwining relation i?(i,i)(5±,i<5±,i = Q±,iQ±,iH ■ But the operator Q ±1 Q±,i is 
a polynomial in H, and we get #(1,1) = H. In general, there holds Huj\ = H. 

The other rule tells that the choice of sequent separations is "commutative". Con- 
struct the Hamiltonian H 1 using the supercharge Then, choosing the second sep- 
aration, we construct H/ 2 ,i) with help of operator Q+,2 which annihilates the states at 
the edges of the second gap of the system H 1 . Hamiltonian H/ 2 u satisfies a relation 
H(2,i)Q±,2Q±,i = Q±,2Q±,\H . The operator Q+,3 = Q±,2Q±,i is of the fourth order and can 
be factorized in different ways. For instance, its alternative factorization is Q +j3 = Q±,iQ±,2, 
that corresponds to the interchanged choices of the separations. Speaking in general, there 
holds H( k j) = H(j >k) . 

Let us denote shortly (j) (or {j,k)) a transformation of Crum of the second (or of the 
fourth) order which annihilates singlet states at the edges of the j-th gap (or of the both 
j-th and A;-th gaps). In this notation, Darboux transformation associated with the ground 
state is represented by (0). Coherently, we denote H^ or Hq^ the Hamiltonians obtained 
by these transformations. Then the rules for sequent use of Darboux-Crum transformations 
can be depicted by the following schemes 



where Id. represents an identity operator. 

In the two-gap case, we can find three new isospectral Hamiltonians in this way. Let 
us present isospectral transformations of the (m = 3, 1 — 0) Lame Hamiltonian with their 
relation to the original system. Seven new isospectral Hamiltonains are found and four self- 
isospectral pairs can be formed. In the following scheme any of the new systems can be 
reached by sequent application of Darboux (0) or second-order Crum's transformation (k) 
on the initial Hamiltonian H . The vertical lines correspond to transformation representing 
the natural separation so that the Hamiltonians related by these lines form self-isospectral 
pairs. 
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As it can be observed in Fig. 2, the nature of the potentials of the obtained systems is 
distinct from the original ones. However, the spectrum of the corresponding Hamiltonians is 
completely identical. The potentials can be tuned with modular parameter k that broadens 
the applicability of these systems. The infinite period limit of the self-isospectral extension 
of the associated Lame systems is discussed in detail in the forthcoming section. 



5 Superextended Poschl- Teller system and infinite pe- 
riod limit of tri-supersymmetry 

In this section we analyze in detail the infinite period limit of the tri-supersymmetric self- 
isospectral extension (4.18)-(4.21) of the associated Lame system. We explain how the struc- 
ture of tri-supersymmetry is modified and study its implications, in particular the restoration 
of the unbroken tri-supersymmetry. 

Infinite period limit corresponds to k — > 1 (k' — > 0). In this limit the associated Lame 
Hamiltonian changes to the energy operator of the Poschl- Teller system. The band structure 
transforms in the following way. The states of the conduction band are transformed into the 
states of the scattering sector of the spectrum, and the singlet edge-state of the conduction 
band is transformed into the first (lowest) singlet state of the continuous spectrum. The 
valence bands shrink, two band-edge states corresponding to the same permitted band con- 
verge smoothly in a unique wave function. This wavefunction can be non-physical in one of 
the limit systems. 



5.1 Self-isospectral supersymmetry in the infinite-period limit 

When k tends to one, the Jacobi elliptic functions cease to be doubly periodic as their 
real period extends to infinity while the complex period takes a finite value 2iK' = 2m\ 
they transform into the hyperbolic functions, dnrr — > sechx, cnx — > sechx, sni ^ 
tanhx. In this limit, the superextended setting described by the two mutually shifted periodic 
Hamiltonians (4.18) acquires the following form 

where the resulting operators represent two systems with the Poschl- Teller potential of dif- 
ferent interaction strengths specified by the integers m and /, 

d 2 d 2 
H rn = -— - C m sech 2 x, H+ = -— - C, sech 2 a;. (5.2) 
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The system keeps the parity- invariance, [R, =0. As we deal with integer values of m 
and /, Hamiltonians (5.2) are reflectionless, i.e the transmission coefficients are equal to one. 

The Hamiltonian H~ {Hf) possesses m + 1 (I + 1) singlet states, m (I) of them are 
bound states, the remaining one corresponds to the lowest state in the scattering sector. 
Hamiltonians (5.2) are almost-isospectral. Their spectra coincide in the continuous part 
E E [0, oo) and just in I + 1 singlet states. The Hamiltonian H~ has additional m — I 
lower energy levels. This indicates what happens with the spectral structure of the extended 
Hamiltonian l-i m ,i when we stretch the real period into the infinity. The 2m + 1 band-edge 
states of H^ nl transform into m + 1 physical states of while in the case of the super- 
partner system 1 only 2/ + 1 band-edge states of highest excitations converge to the 
physical wave functions, the rest is physically unacceptable. Thus the doublets of band-edge 
states and quadruplets of the quasi-periodic states change into m— I singlets and Z + l doublet 
states, and into the quadruplets of the scattering states (see Fig. 3). The presence of the 
singlet states gives a taste of a different nature of the tri-supersymmetry which we discuss 
in what follows. 

The Hamiltonians (5.2) admit the following representation 

H m = -V_ m V m -m 2 , H^-V^-l 2 , (5.3) 
where the definition and the basic properties of the operator V n are 

V n = D + ntanhx, V\ = -£L n , V_ n V n = V n+1 V_ n ^ + (2n + 1). (5.4) 

We shall focus to the properties of tri-supersymmetry and of the supercharges (4.19)-(4.21) 
in particular. Taking the limit k — 1, the non-diagonal integrals (4.20), (4.21) transform as 
follows 

where the non-diagonal components are 

V^V h (5.7) 

X m ,l = • • • ^m-l^m) ^m,l = ^-mV-m+1 ■ ■ ■ ^-i-2^-!-l, (5.8) 

andX^ = (-l) m -'(X m i )t, = (-l) m+,+1 (Y^ > ,) t - The limit does not violate commutation 
relations (tri-supersymmetry is maintained) and we have 

[7~^t3> ^mj] = \)^n,li 3^mJ] = [^mj> ^mj] = 0- (5-9) 

The components of the squares of the non-diagonal supercharges, 

«T) 2 = ( Xm 'f m " * A_ ) , 0® 2 = f Ym { : > 1 - V ) , (5.10) 

U A m,l A m,l / V U r mJ r m,i 
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correspond to the integrals of motion of individual Poschl- Teller subsystems. The results of 
the previous section suggest that these will be proportional to certain spectral type polyno- 
mials in Hamiltonian. This is the case indeed. However, the situation changes significantly 
comparing with the periodic system. With the sequent use of the identity in (5.4), we can 
derive 

m— I — 1 

«T) 2 = I! «J - E mj) = P P X T (KZ)- (5.11) 
3=0 

Here E m j = — (m — j) 2 , j = 0, . . . , m — I — 1, correspond to the m — I singlet states of 
the superextended system (m — I bound-states correspond to the lowest energies of H~ or, 
equivalently to the m — I nonphysical states of H^). The square of the second non-diagonal 
supercharge can be factorized with use of Px T 

m—l 

(y%) 2 = {HZ - E m , m ) n {nZ - E^fpnnZl (5-12) 

j=m-l 

Apart from the roots shared with Px T , there also appear the / + 1 doubly-degenerate energies 
of the superextended system (5.1). The / bound-states energies E m j = — (m — j) 2 , j = 
m — m — 1 , are the double roots, and the energy of the lowest state of the continuous 
spectrum, E m>m = 0, is a simple root. 

Considering the limit case of the diagonal supercharge (4.19), we encounter an interesting 
situation. In [64] it was observed that for a reflectionless Poschl- Teller (PT) system, there 
exists a hidden bosonized supersymmetry. If the system has n bound states (and hence n + 1 
singlet states), there is a parity-odd integral of motion of order 2n + 1, 

V n , (5.13) 

which annihilates all the singlet states and some non-physical states, whose origin was clar- 
ified in [65] from the point of view of the Lame equation and its hidden bosonized super- 
symmetry 9 . The subsystems H~ and have odd-integrals of motion Ai m +\ and A21+1 
of orders 2m + 1 and 21 + 1, respectively. On the other hand, we know that the diagonal 
components of integral Z PT are the parity-odd integrals of order 2m + 1 for each subsys- 
tem. Particularly, would have two parity-odd symmetries of different orders. Let us 
explain the picture and show how the integrals A2 m +i an A21+1 manifest their presence in 
the tri-supersymmetric scheme. 

In the limit Z m ^i — > we can trace the presence of the integrals ^ 2m +i an d A21+1 

k=l ' 

in the diagonal components 

Z rn,l Z m,l = X m,l Y m,l = ^2m+l = V- m V^ m+1 . . . V . . . V m ^{D m , (5.14) 

Z m,l Z m,l = X m,l Y m,l = X m,l X m,l^ 2l + 1 - (5.15) 

9 Earlier, higher order differential operators of this nature were discussed in the context of supersymmetric 
quantum mechanics in [48, 66, 67]. However, their sense and the intimate relation with the algebro-geometric 
potentials were not understood, see also footnote 7. 
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Each of these operators (it is worth to note again that they have the same order) is the integral 
for the corresponding subsystem, and together they annihilate the singlet and doublet states 
of the super-extended system. (>2^J) 2 produces a polynomial of the form 

m— 1 

{Z P Zf = (HZ - E m , m ) J] {HZ - E m ,)\ (5.16) 

j=0 

which can be related with a degenerate spectral hyperelliptic curve of genus m, in contrary to 
the m-gap system (4.18), whose spectral polynomial (3.22) corresponds to a non-degenerate 
hyperelliptic curve of the same genus. It reflects the fact that the band structure disappeared; 
every two band-edge states of the same band transform into a single bound state which ends 
up in the degeneracy in the spectral polynomial. Accordingly, the degeneracy does not 
appear for the lowest state of the continuous spectrum. 

The components of the integral can be rewritten in the following way 

Y m,i = A 2 i+iDi + i . . . V m _{D m = A 2 i+iX~, Y+j = X J h A 2 i+i. (5.17) 
The relation (5.12) can be expressed also as 

(yS) 2 = AlunZ)pnnZ), (s.is) 

where in correspondence with (5.13) A21+1 = T>_{D_i +1 . . . T> . . . T> \-\D x . These formulas 
provide an alternative insight into the kernel of the supercharge 3^T anc ^ commuta- 
tion relation with T~tZ\- The later one can be derived independently just with use of the 
commutation relations [Ti-Z^mJ] = an d W^,A 2 i + ±] = 0. In particular, we have 

Y m,i H m = Ati+iX-jH' = A 2i+1 HfX m j = HfA 2 i +1 X ml = HfY'j (5.19) 

In the infinite-period limit, the tight relation of the non-diagonal integrals (5.5) and (5.6) is 
manifested by means of the parity-odd integral A21+1 of Hf, see Eq. (5.17), which was not 
presented in the periodic case. As a consequence, this integral appears also in the structure 
of the diagonal integral Z*^ see (5.15). 

5.2 Supercharges action and relation of nonphysical with physical 
solutions 

In the case of associated Lame system, the action of the operators l and l was quite 
clear. Each of them annihilated two disjoint subsets of the whole family of 2m + 1 singlet 
states of the Hamiltonian H~ v In the limit case, the situation ceases to be so transparent. 
The systems described by (5.2) differ in the number of the singlet states, the Hamiltonian 
H^Hf) has m + 1 (/ + 1) bound-states. However, the order of the supercharges is not 
affected by the limit. There arises a natural question: what kind of functions is annihilated 
additionally by the nontrivial integrals. We clarify here this intricate situation. 

We introduced the operator V n (see (5.4)) which proved to be useful in factorization of 
both the Hamiltonians (5.3) and the supercharges (5.5), (5.6). We can interpret this operator 
as a Darboux tranformation which satisfies 

V m H~ = H-_{D m , (5.20) 
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where the new Hamiltonian H m _\ has m — 1 bound states. The operator V m annihilates 
the ground state of H~ so that the corresponding energy level is missing in the spectrum 
of H m -\. The Hamiltonians H m and H m -i related by T> m are of the same nature but with 
a shifted parameter. This phenomenon, mediated by Darboux transformation, is called a 
shape invariance. 

We can apply the transformation of Darboux repeatedly, annihilating the lowest bound 
state in each step. This procedure induces the following sequence of Hamiltonians 

H m H m-1 H m-2 •••• ^ H l+1 ^ H r = H h ( 5 - 21 ) 

The resulting operator which relates H~ and can be interpreted as a Crum-Darboux 
transformation of order m — I. This transformation coincides with that produced by the 
operator X~ x 

x m ,iH m — Hi X m l , X m j = T>i + i . . . V m . (5.22) 

Therefore, X^ ml annihilates m — I bound states of H~. In the sense of the superextended 
Hamiltonian 7i m ^, these bound-states are singlets. The energy levels corresponding to the 
states annihilated by X m l are absent in the spectrum of (that makes the systems almost- 

isospectral). Let us denote the corresponding m — I bound-state energies of H~ as E m . We 
denote I + 1 remaining energies of singlet states of H~ as E t . These I + 1 levels are shared 
by the singlet states of Hf. So, the energies of singlet states of H~ are formed by E m and 

The operator X^ n t annihilates a non-physical eigenstate of corresponding to the 
energy E m = 0, Hf-ip m = E m ijj m Acting with the operator 1 on the / + 1 phys- 

ical states (/ bound states and the lowest state of the continuous spectrum) ipi of Hf, 
Hftpi — Eitpi, we get the bound states of H~ corresponding to the energy E t X+^i = ipi, 
H~ipi = Eiipi. The operator Y* t annihilates the physical eigenstates ipi, Y^i^i = 0? an d 
also annihilates the second, non-physical, solution fj m of the Hamiltonian correspond- 
ing to the eigenvalue E m , Y+fi m = 0, Hfrj m = E m fj m , fj m = ip m f x $~ 2 dx. The function 
f\ m is mapped by X^ to a physical state ip m of H~ with energy E m . The same result is 
obtained when we act by on the state ijj m G KerX+,, Y+$ m = ifj m , X+ nl f\ m oc ip m , 

H~i> m = E m ^) m . 

To get an insight into the kernel of the operator Y ml , it is convenient to rewrite this 
operator in the factorized form 

Y-j = V-iV-i+i ...V-i VqVi ■ ■ .Vm-iVm - (5.23) 

annihilates singlets 

The indicated right part of the operator annihilates all the m + 1 singlet states of H~. Due 
to the remaining part, there are additional annihilated / functions <f>j, j = m — /,..., m — 1. 
These additional functions need not to be solutions of the Schrodinger equation corresponding 
to H~. Indeed, due to representation 

m—1 

KKi = ft- - E mjm ) J] (H~ - E md ) 2 P^ T (H-), (5.24) 

j=m-l 
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the functions <f>j satisfy (H m — E m j)(f>j = ipj, (H m — E m j)ipj = 0. Hence, the Hamiltonian H m 
restricted on the kernel of Y ml contains Jordan blocks associated with the energies E t . Due 

to the similar structure of the spectral polynomial of (5.16), the non-physical states 
annihilated by the supercharge of the hidden bosonized supersymmetry are of the same 
nature. 

It is instructive to discuss the case of I = in detail. In this case Hq corresponds to the 
free particle. The commutation relations (5.9) tells that the free particle energy operator is 
related with H~ via the following intertwining relations 

■X-mflHo = H m X^ , Y+ q Hq = H m Y+ fi , (5.25) 

X rn flH m = H^X mQ , Y m0 H m = H^Y m0 . (5.26) 

The first relation of (5.26) is mediated by X~ Q = V\...V m . Keeping in mind (5.22), the 
supercharge X~ annihilates all the bound-states of H~. In the second relation of (5.26), 
the Hamiltonians are intertwined by Y~ = T>oD\...T> m . Apparently, this operator makes 
the same job as X~ but Y~ annihilates additionally the first scattering state of H~. 

The operator X~ transforms the lowest scattering state of H~ into a constant function, 
the scattering state of a free particle corresponding to the lowest energy. This function is 
annihilated by Y* Q while applying X^ we get the initial first scattering state of H~. In 
general, the operators X^ and Y^ transform solutions of Schrodinger equation correspond- 
ing to the free particle Hq into the (formal) eigenstates of H~. They can be employed in 
reconstruction of the scattering states of the Hamiltonian H~ from the plane wave states of 
a free particle, 



i>± = X+ t0 e ±tKX = V_ m V_ m+1 ...V_ 2 V_ 1 e ±lKX 



oc y^ ±iKX = V^ m V_ m+1 . . . V_{D^\ (5.27) 



where ip^ = e ±tKX satisfies 



= Ej±, H~^ = E K ^, E k = n 2 . (5.28) 

Let us summarize the obtained results. Extending the real period of the self-isospectral 
extension of associated Lame Hamiltonian to infinity, the associated superalgebraic structure 
was modified. The squares of the supercharges Z^ T t and turned out to be degenerated 
polynomials in T~t^- This modifies the structure of the underlying superalgebra of the 

bosonic operators in the dependence on a chosen grading operator r*. Comparing with 
the periodic case, the bosonic operators form the algebra u(l) © e(2) © e(2) for the singlet 
energy levels. Recall that the periodic extended tri-supersymmetric system does not have 
singlet states in its spectrum. 

The operators X ml and Y ml annihilate the m — I lowest bound states of the Hamilto- 
nian H~. The eigenvectors of corresponding to these energies cease to be physically 
acceptable. Consequently, isospectrality of the initial system is broken. Speaking in terms 
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of the extended system, the supercharge X£J annihilates the singlet states. The operator 
y^ T annihilates both doublets and singlets which are annihilated by the diagonal operator 
Z^ T as well. From this point of view, the spontaneously (or, dynamically) partially broken 
tri-supersymmetry of the periodic system is recovered in the infinite period limit. 

6 Concluding remarks and outlook 

In the particular case of associated Lame systems, the results of the present paper should be 
understood in a broader context of the existing literature. Dunne and Feinberg considered 
the class I = m—1 of associated Lame Hamiltonians (4.2) as an example of the self-isospectral 
extension provided by Darboux transformation [25]. Khare and Shukhatme [27] found that 
this transformation provides a self-isospectral extension of pure Lame systems just in the 
one-gap case while for the other setting the extension proved to be of a completely different 
nature. On the other hand, Fernandez et al revealed self-isospectrality of two-gap Lame 
Hamiltonian when the second-order transformation was applied [30]. 

In the light of the presented results, we can understand those findings just as pieces of the 
mosaic, which was fully unfolded by the structure of the tri-supersymmetry and especially by 
the self-isospectral supersymmetry of the associated Lame system. In particular, the system 
considered by Dunne and Feinberg is the self-isospectral extension l-L m ,m-\ of the associated 
Lame Hamiltonian, see (4.21). Besides the first-order supercharge A' m m _ 1 , the list of its local 
integrals of motions should be completed by the other non-diagonal supercharge y mjm -i and 
diagonal integral Z m m „i which plays the role of the central charge of the resulting extended 
N = 4 nonlinear supersymmetry. Although both X m ,m-i and Z m ^ m -\ annihilate the doublet 
of ground states, the tri-supersymmetry is spontaneously partially broken since the doublet 
of ground states does not vanish under the action of the supercharge y m ,m-i- This suggests 
that the supersymmetry breaking should be analyzed having in mind the complete set of 
nontrivial local integrals, which are Z, and in the case of the studied general class 
of finite-gap systems. 

On our way to the presented results we left untouched various appealing questions and 
problems. For instance, the self-isospectrality conjecture could be tested on the finite-gap 
systems with missing anti-periodic states. Since these should be prevented from the self- 
isospectral extensions, the structure of the tri-supersymmetry could exhibit peculiarities in 
this case. Besides, the exact proof of the conjecture should be provided. 

The infinite-period limit could be an effective technique in production of the tri- 
supersymmetric systems with non-periodic potentials. In the limit case of the self-isospectral 
extension of the associated Lame Hamiltonian, the isospectrality was broken followed by the 
recovery of the exact tri-supersymmetry. There appears a natural question whether this is 
the common feature or there exist isospectral extensions of non-periodic systems with broken 
tri-supersymmetry. The limit of other isospectral extensions of associated Lame system could 
provide an insight into the general situation. The relation of the tri-supersymmetry and the 
representations of Lie algebras might give an interesting insight into involved physical models 
as well. 

Our construction of the tri-supersymmetric extensions was based on the specific factor- 
ization of the odd-order integral of motion. Relaxing the smoothness of the potential, the 
formal construction should be applicable on the broad family of algebro-geometric poten- 
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tials where the presence of the parity-odd diagonal integral Z is guaranteed. It is worth to 
mention the Treibich-Verdier family of potentials [68] in this context. Besides the associated 
Lame systems, this family contains singular potentials, which could be convenient examples 
to study the tri-supersymmetry in singular systems. 

Regular Crum-Darboux transformations with zero modes in the prohibited bands can 
produce self-isospectral potentials with a generic shift of the coordinate, or superpartners 
with periodicity defects. The particular results of this type were obtained in [28, 29, 44, 45, 
46] with making use of the first- and the second-order transformations applied to one- and 
two-gap Lame equation. It would be interesting to analyse such a class of systems on the 
presence of the tri-supersymmetric structure. 

The revealed super symmetric structure was based on the internal properties of the inte- 
gral of motion Z, related with the KdV hierarchy. This indication of the tri-supersymmetry 
and self-isospectrality in the context of nonlinear integrable systems should be followed 
and analyzed. A special attention should be paid to possible manifestations of the tri- 
supersymmetry in physical systems [34]. 

The work has been partially supported by CONICYT, DICYT (USACH) and by FONDE- 
CYT under grants 1050001 and 3085013. We are grateful to V. Enolskii, V. Spiridonov, A. 
Treibich, R. Weikard and A. Zabrodin for valuable communications. Our special thanks are 
to B. Dubrovin for many detailed explanations on the theory of finite-gap systems. 



Appendix A 

Higher-order differential operators play the key role in the construction of the tri- 
supersymmetry since they mediate intertwining of the super-partner Hamiltonians. As we 
explained in the section on the Crum-Darboux transformations, properties of these operators 
are determining for the physical characteristics of the superpartner systems. We present here 
a short resume of the relevant facts referring for the details to [38], [40]. 

Consider a differential operator of order n which annihilates n functions ipi, i — 1, . . . , n, 

n 

A n = D n + J24( x ) Dn ~ j > 4,^ = 0, i = l,...,ra. (A.l) 

3=1 

Its coefficients are determined by the functions ipi. For instance, the coefficient cf(x) can be 
given in terms of the Wronskian of the n functions ipi, cf(x) = — 4- In W(ipi, . . . , ip n ), where 

W(ipi, . . . , ip n ) = W = det B, Bij = d \ , i,j,= 1, • • • , n. This is in accordance with the 

general form for the coefficients c^(x) = — -j^-, j — 1, • • • ,n, where Wj is the determinant 

of the matrix B modified by replacing the line ^ n--J '\ . . . , ip n n ~^ by ip^i\ . . . , ipn™' 1 ■ In this 
notation, W = W. 

The operator A n can be factorized in terms of the first order differential operators. There 
follow equivalent representations of A n which provide a better insight into the properties of 
the operator, see [38], 

A - ( 1)" Wn D WLl D D^L D ^ (A 2) 
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We can write equivalently 



A n = L n L n _i . . .L 2 Li, Lj=D-aj, aj = ^\n-^-, j = l,...,n. 
The operator can also be expressed as a determinant 



(A.3) 



V>1 


^2 

^2 " 


• < 


1 

D 


^1 


# • 


• ^ 


D n-1 
D n 



(A.4) 



where the multiplicative factor fixes the coefficient of D n to be equal to one. Here, 
the determinant of the operator-valued (n + 1) x (n + 1) matrix is defined as detC = 
Y.aeG n+1 s 9 n ( (T ) C a(i),iCa(2),2 • • • C^(n+i),n+i , where G n+ i is a set of all possible permutations 
of the integers {1, . . . , n + 1}. 

Particularly, when ipi are periodic functions except even number of antiperiodic ones, the 
Wronskian W is periodic. Since the derivatives do not change the period of the functions, 
Wi are periodic as well. The formulas above then justify periodicity of the operator A n . 

Finally, let us make a few comments on the superpartner Hamiltonians H and H inter- 
twined by operator A n (see (3.10)) which annihilates a part of the physical states of H. Let 
the potential of H be smooth and the Wronskian W computed on the kernel of A n be a 
nodeless function. Then the potential of H is smooth as well. The operator A n can be used 
in reconstruction of the eigenstates ip of H corresponding to the eingenstates ip ^ ipi of H 
with the same eigenvalue, Hip = Eip, Hip = Eip, ip = A n ip. These wave functions ip can 
be also represented as 



ip = A n ip = 



W{ip u ...,ip n ,ip) 
W{^...,ip n ) 



(A.5) 



This receipt fails in the reconstruction of the states ipi, which correspond to the same eigen- 
value as ipi, Hipi = Eiipi, Hipi = Eiipi, where ipi is annihilated by A n . These functions ipi, 
annihilated by A^, are given by 



A = 



W(lp!,..,1pi,..,1p r , 



w{ip 1 ,...,ip n ) ' 

where the entry below a symbol is omitted. 



l,..,n, 4^ = 0, 



(A.6) 



Appendix B 

Grading T* = <r 3 

In this case, which corresponds to the usual choice of the grading operator, the non-diagonal 
supercharges Q± are fermionic operators, {Q±, 03} = 0, whereas the diagonal integral Z is a 
bosonic generator. Table 1 represents the explicit identification of the bosonic and fermionic 
generators, and corresponding polynomials appearing in the (anti)commutation relations. 
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Fermionic 
integrals 


Fi = Q- 


F 2 = -Ra 3 Q- 


F 3 = -iRQ_ 


F 4 = ia 3 Q_ 


F 5 = RQ + 


F 6 = -Q+ 


F 7 = ia 3 Q + 


F 8 = -iRa 3 Q+ 


Bosonic 
integrals 


n 


Si = -R 


T* = <J 3 


E 2 = -R<j 3 


B\ = —iRo 3 Z 


B 2 = ~cr 3 Z 


B 3 = -iRZ 


B 4 = -Z 


Polynomials 


P22 = P2B = P-{H) 


P u = P 1B = P + {H) 


Pl2 = l 


Pb = Pz{H) 



Table 1: Integrals of motion and structure polynomials, grading r* = a 3 . 
Grading r* = R 

For the choice = R, the parity-odd diagonal, Z, and non-diagonal, <2_, integrals are 
fermionic supercharges. The non-diagonal parity-even integral Q + is identified as a bosonic 
generator. The identification of all the generators and structure polynomials are given by 
Table 2. 



Fermionic 
integrals 


F 1 = Q- 


F 2 = Ra 3 Q_ 


F 3 = ia 3 Q- 


F A = iRQ_ 


F 5 = iRZ 


Fq = iRa 3 Z 


F 7 = a 3 Z 


F 8 = Z 


Bosonic 
integrals 


n 


Y ^ = R 


Si = a 3 


S 2 = Ra 3 


B 1 = Q+ 


B 2 = icr 3 Q + 


B 3 = RQ + 


B 4 = iRa 3 Q + 


Polynomials 


P22 = P12 = P-{H) 


Pu = Pz{H) 


P2B = l 


Pb = Pis = P+{H) 



Table 2: Integrals of motion and structure polynomials, grading = R. 



Grading T* = Rcr 3 

With this choice of the grading operator, integrals Z and Q + are identified as fermionic 
supercharges, integral <2_ is a bosonic generator. Complete identification of the generators 
and structure polynomials are represented by Table 3. 



Fermionic 
integrals 


Fi — Z 


F 2 = -a 3 Z 


F 3 = -iRZ 


F 4 = iRa 3 Z 


F 5 = RQ + 


F 6 = -Q + 


F 7 = iRa 3 Q + 


F 8 = -ia 3 Q+ 


Bosonic 
integrals 


n 


Si = -R 


S2 = — <J 3 


r, = Ra 3 


B 1 = -ia 3 Q_ 


B 2 = -Ra 3 Q_ 


B 3 = -iRQ_ 


B A = -Q_ 


Polynomials 


P22 = Pz{H) 


Pu = P12 = P+(H) 


Pib = 1 


Pb = P2B = P-(H) 



Table 3: Integrals of motion and structure polynomials, grading = Ra 3 . 

The anti-commutation relations between the fermionic operators are given in Table 4, 
while Table 5 provides the boson-fermion commutation relations. 
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Pi 


F 2 


F 3 


P4 


P> 


P> 


F 7 




P 


P22 


S2-P22 











B A P 12 





PlPl2 


F 2 


S2-P22 


P22 








-B 2 P 12 





-P3P12 





F 3 








P22 


S2-P2 





-P3P2 





P2P2 


F 4 








S2-P22 


P22 


P1P2 





-P4P12 





P 5 





"P 2 Pl2 





SlPl2 


Pu 










F 6 


P4P2 





-B 3 P 12 







P11 








F 7 





-P3P2 





-P4P2 








P11 




F 8 


SlPl2 





5 2 Pl2 











SiPn 


Pll 



Table 4: Fermion-fermion anti-commutation relations. Here the overall multiplicative fac- 
tor 2 is omitted. To get anti-commutator, the corresponding entry should be multiplied by 2, 
for instance, {Pi, Pi} = 2P 22 . 





1 F i 


P 2 


P 3 


P 4 


P 5 


P 6 


P7 


Ps 




-iF 4 


-iF 3 


iP 2 


iPi 


iF 8 


iF 7 


-*F 6 


-*p 5 


£1 


-iF 3 


-iF 4 




?P 2 














£ 2 














iF 7 


tF 8 


-iP 5 


-iF 6 


Pi 





tF 6 P 2B 


-^P 7 P 2B 








-iF 2 P 1B 


*P 3 Pl£ 





p 2 


iF 7 P 2B 








iF^P 2B 





-iF 4 P lB 


-iFiPib 





B 3 


-iF 5 P 2B 








iF % P 2B 










-tF 4 P 1B 


P 4 





-*F 8 P 2B 


-*F 5 P 2B 





tF 3 P lB 








%F 2 P lB 



Table 5: Boson-fermion commutation relations. The overall multiplicative factor 2 is omit- 
ted. To get commutator, the corresponding entry should be multiplied by 2, for instance, 

[r„p] = -2*f 4 . 



Appendix C 

In the treatment of the section 4.1, we left untouched the system described by Lame associ- 
ated Hamiltonian with m — l. Due to an identity dn (x + K) = k'/dnx, in contrary to the 
other members of the family, its period is K. This fact explains why the algebraic methods 
applied to other members of the family in this case say that the dimension of sl(2, R) repre- 
sentation realized on antiperiodic in the period 2P singlet states is equal to m — I — 0. This 
is just because singlet states with such a period does not exist. The place of this system 
in the mosaic of the tri-supersymmetric self-isospectral systems is clarified by its intimate 
relation to pure Lame system, mediated by Landen transformation [69, 70]. 
Landen's transformation of the elliptic functions can be written as 

sn (-, k) cn (-, k) 1 — asn 2 (-,/t) 
sn (x,k) = a a ' \ a — '-, cn (x, k) — — - a , 

, x «' + (2 - a) cn 2 (^, k) , x 

where a = k 2 = jp^p. , k = j^t, a — 1 + k'. To avoid confusions, let us denote explicitly 
the dependence of the complete elliptic integral K on the modular parameter such that we 
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will write K{k) or K(k). Since K(k) = (1 + k)K{k), Landen's transformation divides in 
two the period 2K(k) of the elliptic functions in the sense that the period of the resulting 
expression is K(k). 

Using the identities (C.I), we can rewrite the Lame Hamiltonian in terms of the el- 
liptic functions of a new variable y — ^ and the modular parameter k, H~ {x,k) = 
^2 [H mjm (y, re)] + const. The displacement K(k) of the pure Lame tri-supersymmetric partner 
changes to in the case m — I. It is in accordance with our result on the general case 
m 7^ I where the superpartner potential was displaced in the half of the real period. Thus, 
we obtain finally the relation 



where c is a constant term. 

It suggests directly the form of the tri-supersymmetry in the special case of m — I Lame 
associated systems; all the operators commuting with H m fi commute with H mjTn as well. 
To get their explicit form for the systems described by 7i m ,m we just have to rescale the 
variable and apply the identities (C.l) in the formulas (4.21) for X m> o(x, k), y m fi{x,k) and 
Z m fi{x, k). Then we can write immediately 

X m ,m(y, re) = Xmfl {ay, k{n)) = X m>0 ^(1 + n')y, j - ^ , 

y m ,m{y,K) = y mfi {ay, k{n)) = y mfi ^(1 + K r )y, | - ^ , 

Z m , m {y, re) = Z mfi {ay, fc(/c)) = Z mfi ^(1 + n')y, j - ^ . (C.3) 

Naturally, the algebraic relations between the operators remain unchanged. Thus, the self- 
isospectral supersymmetry and associated superalgebra are recovered for m — I case. 
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Figure 2: The Hamiltonians H^^) = —D 2 + V^-m — 12 with the plotted potentials are 
spectrally identical (see the last diagram in sec. 4). The potentials on the left and right are 
mutually shifted in half of the period. The modular parameter is set k = 0.99. 




m = 5, / = 3 



Figure 3: Potentials on the left correspond to H~ nl (dashed thick line) and H^ tl (solid thick 
line); k 2 = 0.99. On the right, the limit k — > 1 of these potential functions is shown. The 
solid thin lines represent the shared bound-states and the lowest scattering state. Dashed 
thin lines represent m — I = 2 singlet states 
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